Analogy of the slow dynamics between the supercooled liquid and the supercooled 
plastic crystal states of difluorotetrachloroethane 
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Slow dynamics of difluorotetrachloroethane in both supercooled plastic crystal and supercooled 
liquid states have been investigated from Molecular Dynamics simulations. The temperature and 
wave- vector dependence of collective dynamics in both states are probed using coherent dynamical 
scattering functions S(Q,t). Our results confirm the strong analogy between molecular liquids and 
plastic crystals for which a-relaxation times and non-ergodicity parameters are controlled by the 
non trivial static correlations S(Q) as predicted by the Mode Coupling Theory. The use of infinitely 
thin needles distributed on a lattice as model of plastic crystals is discussed. 
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It is well known that very different materials such as: 
silica, low-molecular-weight liquids (carbohydrates, alco- 
hols), polymers or even proteins are able to exhibit a 
very intriguing feature called glass transition 0. This 
latter is characterized by an extraordinary decrease of 
several orders of magnitude of the mobility in a narrow 
temperature range. Translational, orientational or even 
internal degrees of freedom (TDOF, ODOF or IDOF) 
can be at the origin of the mobility of a specific sub- 
stance and approaching the glass transition transition 
temperature T g , all of them are inextricably coupled. 
Most of the investigations performed on glassy materi- 
als considered molecular liquids and focused on TDOF 
so far. The precise role of the ODOF remains par- 
ticularly unclear as demonstrated by the so-called low- 
temperature translation-rotation paradox Q. Plastic 
crystals i.e molecular crystals composed of orientation- 
ally disordered molecules offer an interesting solution in 
order to mainly investigate ODOF during the freezing 
process. Indeed, some of the plastic crystals called glassy 
crystals 0] such as cyanoadamantane ethanol j]| or 
orthocarborane Q can be considered are true rotational 
analogs of canonical liquid glassformers since they show 
a step of the specific heat at T g or a non-Arrhenius be- 
haviour of the rotational relaxation times. Only a few 
substances show the extraordinary property to exhibit 
a glass transition from both the plastic and the liquid 
phase: ethanol cyclohexene @, PMS or 2- 

bromothiophene [l0(. An investigation of such materials 
is particularly interesting to understand the fundamental 
microscopic mechanisms of the glass transition and the 
precise interplay between the different degrees of free- 
dom. 

In the high temperatures pico-nanosecond regime (ps- 
ns), a fundamental question concerns the onset of the 
precursor cooperative mechanisms which lead to the glass 
transition at T g . In the last two decades, much at- 
tention has been devoted to the mode coupling theory 
(MCT) 0. So far, it is the only theory which provides 
a microscopic description of supercooled atomic liquids. 



The intrinsic basis of MCT states that the behaviour of 
any time-dependent correlators describing the dynamics 
of the system is only controlled by its static density cor- 
relator i.e S(Q) and predicts the existence of a critical 
temperature T c . Scaling properties predicted by MCT 
approaching T c have been successfully validated in nu- 
merous experiments and molecular dynamics (MD) sim- 
ulations |l2| of systems whose dynamics are controlled by 
TDOF. Recently, some extensions of MCT called molec- 
ular mode coupling theory (MMCT) 0] have been pro- 
posed to take ODOF into account: one diatomic probe 
molecule in an atomic liquid, liquids made of linear 
molecules or water to cite only a few. The authors have 
particularly shown that some of the basic predictions of 
MCT still hold owing TDOF/ODOF coupling. As re- 
vealed by a recent MD investigations [ljj performed on 
orthotcrphenyl (OTP), coupling of the rotational dynam- 
ics to the center-of-mass motion can be very complex. No 
microscopic theory has been developed for plastic crystals 
so far. However, in [3. Il4|. we have particularly shown 
from NMR and Raman experiments, and MD computer 
simulations that some predictions of the idealized version 
of the MCT (critical temperature T c and time scaling 
laws) were able to describe relatively well rotational dy- 
namics of different plastic crystals. This intriguing previ- 
ous result call for new investigations to clarify the similar- 
ity between slow dynamics behaviour of plastic crystals 
and molecular liquids. 

In this Letter, we present results of a Molecular Dy- 
namics (MD) comparative numerical study of the su- 
percooled plastic and the supercooled liquid phases of 
difluorotetrachloroethane (DFTCE). It should be noted 
that the supercooled state of liquid DFTCE is ob- 
tained owing the hyperquenching rate of the MD sim- 
ulation. This compound is composed of simple molecules 
CFCI2 — CFCI2 close to dumbbells extensively used in 
MD calculations as prototype of molecular glass-formers 
liquids DFTCE plastic crystal has been experi- 

mentally widely studied and exhibits a glass transition 
of the overall rotation of the molecules at T g = 86 K . 
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Changes concerning the nature of dynamics in this sys- 
tem has been reported from NMR experiments, Brillouin 
and dielectric spectroscopy where it was suggested that 
the freezing process could be described on the basis of 
MCT Hi . 
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FIG. 1: Coherent intermediate scattering function Sd{Q,t) at 
Q — 1.54 A obtained for the plastic phase from T = 130 to 
260 K and the liquid phase from T = 150 to 300 K. 

MD simulations were performed for an orientationally 
disordered crystal sample of DFTCE at 14 different tem- 
peratures from T = 130 to 260 K in steps of 10 K and 
a supercooled liquid sample at 13 different temperatures 
from T = 150 to the melting temperature T m ~ 300 K in 
steps of 10 K. The same model was used for both phases 
and it is completely described in [rR Hsf . 

The intermediate scattering function as it can be clas- 
sically obtained from coherent neutron scattering experi- 
ments will be mainly considered: S(Q,t) — (Pq(*)Pq(0)) 
where Pq{t) is the time-dependent density correlator: 

Pq{t) = b a exp[iQ.r a (t)] where the sum is over all 
the atoms a of the system. b a and f a are the coher- 
ent scattering length and the position of the a atom 
respectively. An average over isotropically distributed 
Q-vectors having the same modulus Q is performed in 
order to obtain S(Q,t) for a polycrystallinc sample. In 
general, S(Q,t) can be expressed as 0: S(Q,t) — 
S C (Q) + Sd(Q,t) where S C (Q) = \{pq~}\ 2 is the coherent 
elastic scattering and Sd(Q,t) = {SpQ(t)5p^(Q)) where 
Sp^(t) — Pq{t) — (Pq) is the fluctuation of the time- 



dependent density operator. For plastic crystals, S(Q,i) 
is identical to Sd{Q, t) except for Q-vectors corresponding 
to Bragg peaks for which the long time limit of S(Q,t) 
reaches the non-zero value S C (Q). This behaviour is 
found at all temperatures and it is associated to the crys- 
talline order of the TDOF and not the freezing process 
of the ODOF. Therefore, in the following, S d (Q,t) will 
be preferred to S(Q,t) in order to be able to compare 
directly molecular liquids and plastic crystal behaviours. 
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FIG. 2: Temperature dependence of the relaxation tq at Q = 
Qmax — 1-27 A for both the plastic and the liquid phase. 
Inverse of the diffusion constant D is also indicated. Solid 
lines show VFT fits (see text). Relaxation times tq at Q = 
0.9, 1.37 and 1.7 A -1 and diffusion constant D (10~ 9 .m 2 s _1 ) 
as function of the rescaled temperature T — T c are displayed 
in inset. Solid lines show fits using power law dependence 
(T - Tc) 1 . (Tj ~ 139 K, 7* ~ 2.55) and (T c p ~ 129 K, j p ~ 
2.45) are found for the supercooled liquid and the supercooled 
plastic crystal respectively. Relaxation times are shifted in 
order to demonstrate the overlap and the generic power law 
dependence. 

Upon cooling, as shown in Fig. ^ both plastic and 
liquid Sd(Q, i) display a two-ste p de cay as classically ob- 
served in supercooled systems [l5j. The long time a- 
relaxation of both plastic and liquid Sd{Q,t) is sepa- 
rated from short time regime (/?) by a plateau-like re- 
gion. This latter can be associated to the so-called cage 
effect corresponding to the trapping of molecules. The 
nature of this cage is of translational origin for liquids 
and of rotational origin (steric hindrance) for plastic crys- 
tals The height of this plateau is usually called the 
non-ergodicity parameter and will be noted fq in the fol- 
lowing. In order to probe the long time regime, we also 
defined the characteristic a-relaxation time tq of the ro- 
tational dynamics as the time it takes for Sd{Q, t) to de- 
cay from 1 to 0.1. In Fig. [3 we plot the relaxation time 
tq at Q = Q m ax corresponding the first diffraction peak 
(see Fig.|3:). In addition, for liquid DFTCE, the transla- 
tional diffusion constant D has also been calculated from 
the mean square displacement (r 2 (t)). tq characteristic 
times and D are found to exhibit a very similar non- 
Arrhenius behavior. In order to quantify this similarity, 
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we used a Vogel-Fulcher-Tammann (VFT) law written 
as t = ^ exp[l/KvFT-(T/TvFT — 1)] where Tvft is 
the temperature of apparent divergence of r and Kypx 
measures the kinetic fragility. We particularly obtained 
K VFT ^ 0.18, 0.16 and 0.20 and T VFT ^ 92.5, 90.4 
and 86.9 K from the inverse of the diffusion coefficient, 
tq of the liquid and tq of the plastic crystal, respec- 
tively. It should be noted that results obtained from MD 
calculations are in fair agreement with the parameters 
KypT = 0.13 and Tvft — 70 K found from dieletric 
experiments of plastic DFTCE [lif . 

In |l4| . we suggested that some common microscopic 
mechanisms, relatively well described by MCT, are in- 
volved in both orientationally disordered crystals and 
molecular liquids. Therefore, in the following, we dis- 
cuss in the MCT framework the wave-vector and tem- 
perature dependence of dynamics for both plastic and 
liquid DFTCE. MCT predictions on the scaling features 
TQ (T) = A(Q).(T - T c )-T in the long time decay have 
been carefully checked using a procedure completely de- 
scribed in [FJj and results are displayed in Fig. For 
plastic crystal DFTCE it was possible to extract the crit- 
ical temperature T£ ~ 129 K and scaling exponent j p ~ 
2.45 from independent fitting analysis performed in the a 
regime at different Q. Using the same method, T\ ~ 139 
K and j l ~ 2.55 were found for the supercooled liquid 
DFTCE. This good agreement with the T-dependence 
MCT predictions is also confirmed by other tests in the 
[3 regime where data were analyzed using a von Schwei- 
dler law including the second order correction which will 
be described elsewhere. Q-dependence MCT predictions 
are presented in the following. A fundamental property 
of MCT or MMCT stresses that dynamics are driven by 
the static density correlator S(Q). This has been con- 
firmed in several numerical and experimental studies for 
simple ,20] or molecular liquids |2j. In order to check 
this behaviour in plastic crystals, we calculated the static 
density correlator S(Q) = S(Q, t = Q) at different tem- 
peratures and wave vectors as well as the non-ergodicity 
parameter /q, the total prefactors h}{Q) and h 2 (Q) and 
the relaxation time tq. A fitting procedure using a von 
Schweidler law as derived from MCT including second or- 
der correction was employed |l4j |. Similarly to the inter- 
mediate scattering function, the static correlator can be 
decomposed into two components: S(Q) = S c (Q)+Sd(Q) 
where S d (Q) = (\6p \ 2 ) = (\ P(3 \ 2 ) - \(p )\ 2 is the diffuse 
scattering. Owing the high rotational disorder, plastics 
crystals yield only a few diffraction peaks but exhibit 
intense and highly structured diffusive scattering which 
provide useful information to characterize rotational mo- 
tions and most probable orientations 0|. Moreover, 
remarkable similarities of the structure factor between 
plastic and liquid phases of different systems have been 
particularly observed 0, H3|- Both S C (Q) and Sd{Q) 
are displayed in Fig. [3: as function of the wave vector Q. 
At the largest Q-vectors, both plastic and liquid Sd(Q) 
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FIG. 3: Q-dependence of different parameters for the plastic 
(full circles) and the supercooled liquid (circles) phase at the 
temperature T = 160 K. a) The a-relaxation times tq /TQ max 
where TQ max is about 345 ps in the plastic phase and 4122 ps 
in the liquid phase at this temperature, b) The non-ergodicity 
parameter fg obtained from a fitting procedure using a von 
Schweidler law as derived from MCT including second order 
corrections, c) The diffuse scattering Sd(Q)- First Bragg 
peaks are also indicated with vertical arrows. 



merge since the intramolecular structure is probed by 
the wave vector Q. tq and /q are also shown in Fig. [3^ 
and b respectively. Both plastic and liquid DFTCE Q- 
dependence clearly exhibit very similar features. In the 
investigated Q-range [0.3 — 5.5] A -1 , S C (Q) shows a few 
Bragg peaks expected for a bec crystalline structure with 
a cell parameter a = 6.82 A -1 (experimentally, a = 6.98 
A -1 ). For Sd(Q), four broad bumps are localized at 1.27, 
2.30, 3.34 and 4.85 A" 1 both in the supercooled plastic 
and liquid phases. The first sharp diffraction peak of the 
liquid particularly corresponds to the first and most in- 
tense Bragg peak of the plastic phase. Clearly, we see 
that /q and tq obtained for plastic or liquid phases are 
mainly correlated with Sd(Q) as predicted by MCT. The 
total prefactors not displayed in Fig. [21 exhibit the same 
behaviour. The non-ergodicity parameter or the total 
pre-factor obtained from S(Q) of the plastic phase also 
shows a modulation with Bragg peaks of the elastic scat- 
tering S C (Q) (not shown in the present Letter). It should 
be mentioned that colliding infinitely thin needles dis- 
tributed on a lattice are often considered as prototype 
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of glassy crystals [i^, 0] • Two-step relaxation processes 
are observed and an ideal glass transition can be obtained 
when the length of the rods (the analog of the temper- 
ature for these systems) reaches a critical value. If the 
thickness of the rods are chosen infinitely thin, all static 
correlations vanish and thus, following MCT predictions, 
no glass transition is expected |U . The present work 
clearly illustrates the strong difference between real plas- 
tic crystals and the colliding needles on the origin of their 
respective ideal glass transition. 

In this Letter, collective dynamics of both supercooled 
plastic crystal and supercooled liquids states of DFTCE 
have been investigated using MD simulations. From the 
calculation of the total intermediate scattering functions 
S(Q, t) at different wave vectors Q and temperatures T, 
our results confirm the strong analogy between two sys- 
tems made of the same molecules but for which dynamics 
are controlled by different degrees of freedom. The au- 
thors hope that the present investigation could contribute 
to develop a theoretical framework for orientationally dis- 
ordered crystal based on MCT 
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